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During  the  year  prior  to  the  contract  with  ONR,  the 
principal  investigator  had  a  similar  contract.  During 
that  year  the  principal  investigator  had  developed  the 
theory  and  algorithms  of  the  prediction  of  the  scattering 
of  light  by  sound  and  the  inversion  procedure  of 
investigating  the  sound  fields  from  the  scattered  optical 
data.  For  the  purpose  of  this  report,  these  theories 
will  he  known  as  the  direct  theory  and  the  inversion 
theory . 

During  the  Spring  of  1980,  the  principal  investigator 
trained  two  students,  Mr.  Charles  Fray  and  Mr.  John  R.  Laflin 
in  the  aspects  of  acousto-optic  interaction. 

The  objective  of  this  contract  was  to  allow  the 
principal  investigator  and  his  graduate  students  to  work 
with  colleagues  in  the  Physical  Acoustics  Branch  of  the 
Naval  Research  Laboratory,  Washington,  DC  to  implement 
algorithms  on  their  computer  data  acquisition  system. 

During  the  stay  at  NRL ,  Mr..  Fray,  in  conjuction 
with  others,  developed  a  modelling  system  which  would 
precidt  the  colored  schlieren  patterns  of  ultrasonic 
fields.  The  output  of  this  model  was  by  colored 
television  display  of  computed  values.  The  work  of 
Mr.  Fray  will  be  continuted  by  those  at  NRL. 
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Mr.  Laflin  pursued  the  inversion  problem.-  That  is, 
he  developed  a  computer  based  experimental  system  to 
acquire  acous'.  o-opt ic  data  and  process  it  to  reveal  the 
complicated  near  field  of  an  ultrasonic  transducer. 

The  principal  investigator  directed  the  students, 
collaborated  on  a  new  theory  for  tomographic  processing 
acousto-optic  data,  and  generally  supported  the  Physical 
Acoustic  Branch  with  theory  and  concepts  in  acousto-optics 
and  scattering  of  sound. 

During  the  extension  period  of  the  contract  from 
September  1980  to  May  1981,  substantial  progress  has 
been  made  in  the  graphic  routines  associated  with  Doth 
the  direct  theory  and  inversion  theory. 

Attachment  A  is  a  preprint  of  a  manuscript  resulting 
from  the  tomographic  work  which  was  presented  f:o  the 
Acoustic  Imaging  conferences,  Monterey,  CA,  Spring  1981. 

At tac hemn t  B  and  C  are  abstracts  of  papers  presented  at 
the  Fall  1980  meeting  of  the  Acoustical  Society  of 
America,  Los  Angelas,  CA. 
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Ml  ABSTRACT  !  ! 

•ss.  j  ! 

30.  The  principles  of  computerized  transverse  tomography  can  be  ! 

JL  applied  to  the  acousto-optic  reconstruction  of  the  local  sound  ■ 
.i?.  pressure  of  an  ultrasonic  field.  |  For  sufficiently  narrow  beams  ; 
.12.  of  ultrasound  in  the  low  megahertz  region,  the  total,  optical  • 
3ji_  phase  retardation  of  an  interrogating  light  beam  can  be  consi-  I 
12.  dered  as  a  projection  of  the  sound  field  pressure.  (Fourier  j 
techniques  for  the  numerical  reconstruction  of  the  pressure  field  • 
3tL  yield  as  intermediate  steps  a  Fourier  domain  associated  with  the  j 
iSL  angular  spectrum  of  plane  waves  comprising  the  sound  field.)  Con-  ' 
n  Isequently  the  sound  field  can  be  reconstructed  in  other  regions  j 
t£L  than  the  plane  of  interrogation,  j  In  this  work  we  discuss  two  al~  ; 
Tl-  |ternative  methods  for  acquiring  data.  One  method  builds  the  • 
jFourier  domain  along  radial  spolces  which  is  inconvenient  for  j 
;numerical  processing  by  DFFT  algorithms.  The  other  procedure  bu~  ; 
lilds  the  Fourier  domain  in  a  nearly  rectangular  format  compatible  ; 
jwith  two-dimensional  DFFT  algorithms.  With  this  latter  method, 
jit  is  possible  to  evaluate  the  pressure,  along  a  line  with  limited  : 
data  and  a  one-dimensional  DFFT.  i  i 

H\  :  i 

ft-  :*  Work  supported  in  part  by  Physical  Acoustics  Center  Program  at 
S&-  !  the  Naval  Research  Laboratory,  Code  5130,  Washington,  D.C. 
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INTRODUCTION 


5^  A  sound  field  of  low  ultrascjnic  power,  low  ultrasonic  fre-  t 

6_  quency,  and  narrow  bean  width  behaves  as  an  optical  phase  grat-  | 
ing.  Collimated  light  passing  through  such  a  sound  field  experi-  j 
8.  ences  an  optical  phase  retardation  proportional  to  the  local  f 
o_  sound  pressure,  integrated  over  the  light  path*  This  constitutes  j 
,0_  j  a  "projection”  of  the  sound  field  and  numerical  methods  of  compu-  j 

i  f erizedl.tr a'nsver'sel.  tomography  can  be _ applied  to _ gs  t iraa te  the  ? 

.2.  local  sound  pressure. _ !_ _ } 

Numerical  techniques  using  Fourier  transforms  are  useful  in  j 
o_  pressure  field  evaluation  since  an  intermediate  step  yields  the 

.o_  Fourier  domairriassociate'd  with  the _ gngular _ spectr.um__0.£ _ pXa.n.e_! 

LJ_  v;aves  comprising  the  sound  field.  Moreover,  hv  modification  of  j 

1 3  the  phase  tefmsO&'f 1  each~ nlanfeP-wave  of  the  angular  spectrum.  it_is  J 

19.  possible  to _ cons truct  the  sourid  field  at  d 3Xf.p.r.en.tL.p lao.e s  .  In  i 

?0  other  words-.  wanAestfina te  of  most  iof  the  sound  field  can  be  com-  ! 
31  puted  from  a  set  of  acousto-joptic  data  taken  over  a  single  • 

32_  transverse,  .plane . _ Ihe_.totaJL_£ieid  >..h.Q.yevexv_ca.nnQt_.b.e..coja£?.truc.tL:-J 

?.3_  ed  everywhere  sinceo.  evahe s'ceritr--waves^-riear^-!l:he)  sound  sx>.uxce_.are_' 
2b_  not  accounted  for.  This  total  field  concept  is  valid  when  the 
25  sound  field  can  be  described  by  'the  Helmholtz  equation,  thus  el~ 
26.  irainating  application  to  non-linear  or  highly  attenuated  sound 
LX.  fields.  | 

| 

L?.  In  a  series  of  papersl-^  directed  toward  transducer  cali- 

30  bration  Cook  and  Berlinghieri  Have  described  one  method  of  data  J 

31  collection  which  we  will  call  Metjhod  A.  Acousto-optic  data  is  j 

32  collected  at  the  terminus  of  tbje  light  paths  as  shown  in  Figure  ; 

33.  1.  These  light  paths  are  parallel  to  each  other  and  are  in  a  j 

S!l  plane  parallel  to  the  surface  of  !the  transducer*  Sufficient  data  J 
3J?_  can  be  acquired  from  interrogation  of  the  field  in  one  direction  I 

3?_  if  the  field  is  symmetric.  If  the  field  is  not  symmetric,  Method 

31  A  involves  rotation  of  the  transducer  about  an  axis  normal  to  the  j 
3.9_  transducer  surface,  such  as  CC* .  ! 

j 

kil  Here,  we  present  an  alternative  method  of  data  collection  j 

1  which  we  will  call  Method  B.  In  Method  B  the  transducer  is  ro-  j 
ill.  tated  about  a  line  AA’  parallel  to  the  transducer  surface  and  lo-  , 
.V.:!  cated  in  the  plane  of  the  light  paths.  The  line  AA’  is  also  per-  j 
J--1-  pendicular  to  the  light  paths.  !  * 

!  i 

r  ■  1 

j  We  will  demonstrate  how  both  methods  allow  the  generation  of  '■ 
|data  in  the  angular  spectrum  (plane  wave  decomposition)  domain  . 
jwith  both  phase  and  amplitude  information  to  allow  evaluation  of  . 
ithe  pressure  field  at  the  plane  of  measurement  using  Fourier  ! 
(transforms. 
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_ Data  ac^€ared^usxn^MethodyA'buird's^:tKdL aSigiilar spectrum  do- 

main  in  a  polar  format  through  a  series  of  one-dimensional  Four¬ 
ier  transforms.  The  pressure  field  at  the  plane  of  measurement 
can  then  be  obtained  by  a  two-  dimensional  inverse  Fourier  trans¬ 
form.  DFFT  algorithms,  however,  require  data  to  be  in  a  rectan¬ 
gular  format.  Two  alternatives  are  to  interpolate  the  rectangu¬ 
lar  data  from  the  polar  data  or  to  perform  a  Hankel 


lawayuroa-  'uieaou.Kxn so interpolation v  o  jura. _ . ii'Ss y i  ut_y2.Ji-.s_ 

[there.  On  the  other  hand.  development  of  efficient  algorithms 
[fbifCHgnkelStfanifBrtas  is  now  an  active  area  of  research. _ 


|  He thod  B  exhibits  three  attractive  features. _ The _ first  is 

that  through.'.T'aHserieslof  one-di’mensional  transforms,  the  angular . 

spectrum  domain  can  be  built  in  a1  format  which  closely _ approxi- 

raa t es  a  rectj'anSular'Vfasterf.l!)IOlThe  second  feature  is  tha t:  ._t:he„dat.a_ 
col lected  by  this  method  lies  midway  between  the  angular. _$peAfc£Uia_ 
domain  and'OIfhel'XtanS-space  domain.  Evaluation  of  the  pressure, 
field  at  the  plane  of  measurement;,  therefore,  requires  only  a  . 
J5jg.rAe5__ojfT(_ji.nve.rse>_o.nerdimei\s.i.clnal.-tr.anaforms.^_._  A._thircLf  eatur.e_ 
is  that  acoustic  pr^ssureican0'6e7cbmpu1ied3^j.'dn^.a^llne_t.r.ansyer.se- 
to  the  direction  of  sound  propagation  by  a  single  inverse 
one-dimensional  transform.  \ 


23.  2.0  THEORY  OF  METHOD  A  i 

22.  | 

30  Consider  a  harmonic  sound  field  being  produced  by  a  planar 
31_  transducer.  Let  the  pressure  at  [a  distance  z  from  the  transducer 
32  be  expressed  as  } 


p(x,y,z,t)  *»  p(x,y,z)exp(-ja)t) 


In  the  following  discussion  the  time  variance  will  be  dropped  for 

convenience.  ! 

i 

I 

l 

Line  integrals  across  the  pressure'  field  at  z+z0  can  be 
written  ! 


p(x,z0)  =  J p(x,y,z0)dy 


Jwhere  the  limits  of  this  integral'  and  others  are  taken  from  minus 
'infinity  to  plus  infinity.  p(x,z0)  can  be  seen  as  a  "projection" 
!of  the  pressure  field  p(x,y,z0).  ! 
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;  _ In  the  !aefeTgrL.ofiu  tHis'Iexperfraant  17,ipCk1rz7>~).'"is\\rdbtained___from  j 

!  measurement  of  the  Raman-Nath' parameter  V  defined  as  | 


V(x,z0)  =  [2Tfk/x3  p(x,z0) 


10  j  where  A  is  the  optical  wavelength  in  vacuum  and  k  is.  the _ medium.  j  10 

11_  ' jYL&zo-frp  tlct coefficient  which  relates  the  index  of  ref r&c .tJkOJX. JL.1 


changes  in  acoustic  pressure.  This  parameter  V  is  a _ measure  of 

£heO optical  bhase-Lxe tarda tion  induced  by  the  sound  field.  It  is 
a  common  paramter  used  in  most  theories  and  can  be  inferred  from 

acous  t  o-o  ptic  measurements. _ 7hLs_J?at£me_teX » ._jn_Qnr_j:as.e.,.._i 

b  e  mea  sur edi'ds^i^haio.c*L_Ya  rioyjs.-te.ch'oiq uea  -S-&X- _ ac  <jui  ring _ the.. 

necessary  pMgg-^d„a,mplitude_.lnfQrm.atlo.n_Qaji_be_XQun<l.ln_.t.he-llrL. 
terature^~QWW^WWPWW>PMl  ! _ _ 

_ I _ _ _ • _ 

We  wil.l'I'shbxD'the^'  relation .  b e_t_^g§ n_t;h e_px£>4 ec t.e.d .  _p r.es s ,ur.e._a nd  j 
the  Fourier  domain  assuming  p(x^zD)  to  be  a  measurable  quantity.  * 
.Sy.b.s.t.i.tLu.tiQn-pJL.a._tw.oa^limensional  .trans  fora  .expression _ into. _ the_r 

illtegral._oJJgqmt^hlIC^)k*giv.e.^L^iiBHH^i^il. - - - J 


p(x,z0)  -  ///p(kx,ky;Zo)expf jClCyXHcyyiJdydlCj-dky  (4) 

! 

i 

» 

where  kx  and  ky  are  components  oi;  the  acoustic  wave  vector  k  . 
The  integration  of  the  y-variable  can  be  completed  yielding  the 
Dirac-6  function  2ir5  (k).  The  j  sifting  properties  of  the 
6  function  upon  integration  over  ;k  produce  the  desired  result 

I 

I 

I 

I 

I 

.  p(x,z0)  =  l/2ir Jp(kx>0;zo)exp(jkxx)dkx  .  (5) 


This  result,  sometimes  referred  ’to  as  the  "Fourier  J 
projection-slice  theorem,”  states  that  the  Fourier  transform  of  a  j 
projection  is  a  slice  of  the  Fourier  transform  of  the  projected  j 
function.  In  other  words,  the  one-dimensional  transform  of  j 
p(x,z0)  produces  a  single  line  in  the  Fourier  domain.  This  line  j 
lies  perpendicular  to  the  direction  of  the  light  paths,  that  is  • 
the  line  of  interrogation.  Consequently,  if  one  rotates  the  j 
transducer  as  specified  for  Method  A  in  equal  angular  increments  ; 
j(essentially  around  the  sound  field  axis),  the  projections  obta-  J 
jined  are  related  to  values  injthe  Fourier  domain  located  along: 
jradial  lines.  In  other  words,  taking  a  discrete  one-dimensional  i 
'transform  of  the  projections  p  (x,z0)  yields  values  shown  as  cir- • 
jcles  in  Figure  2.  This  result  is  not  restricted  to  sound  fields  J 
jand— is-a  -genei'al-  result  of  -  projection-theory. - - - 1 
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_ The  ^revIouslyT  stajedridif f ictil £$jyj\3Ln”  £hJ£jAyfe“ts ion  o£_data_ 

in  the  Fourier  domain  to  the  jtime-space  domain  from  the  polar 
format.  In  addition  to  interpolation  to  a  rectangular  format, 
Mersereau  and  Oppenheira  suggest  other  schemes  to  circumvent  this 
problem. ^  | 


3.0  THEORY  OF  METHOD  B _ ! _ ! _ '  . 

FIRST  0?  TITLE _ 4_ _ _ _ 

_ Consider  the  origin  of  the  Cxry,z)  coordinate  axes  to  be  the,. 

Eoifit'D  o£I?SinEe'rsect:ion  of  the  acoustic  axis_  and_.the_.ix|.ane__Qf_tliff_. 
light  slice.  Again,  z  is  the  acoustic  axis.  The  light  path 
slice  will^piL_1b.e„st_aii_a.n&le_l_Jiith_the.  y-axis* — Jfe._c ajxjyxi t.e_. 
the  projectA9^S?S£heipr^auTta_fi;ald_vstIJtii_th-e.„nat.a.ti.Qii_ch.al»g.e.!i_to„ 
account  for  the  angle  as  _ i _ 


yheicp_dy^5_alQng:\yie:aight_^aths_tatlaTt!:Angl&;<i.Lwlth„the_yriaxls.._ 

I 

The  Fourier  domain  pressure  jean  be  expressed  as 

I 

t 

j 

•  p(kx>ky,z)  ~  Po(,kx>ky)exPOkz55) 

! 

i 

i 

where  pQ(kx,ky)  is  the  Fourier  component  at  z  ■  0  and  fc2  is  the 
z-component  or  the  acoustic  wave  [vector.  He  again  substitute  the 
Fourier  description  of  the  pressure  field  in  Equation  (6)  and  in¬ 
corporate  Equation  (7)  to  give  j 

I 

*  I 

J  . 

p(x,<|>,z)  “ ///P0(kx,ky)exp[  j(kxJs+kyyfkj5z))dy'dkxdky  (8) 

I 

I 

I 

I 

The  light  path  slice  and  the  x-axis  defines  a  nev;  coordinate 
system  (x, y’jZ*).  This  is  related  to  the  (x,y,z)  coordinate  sys¬ 
tem  by  the  transformation 


2 2 

26 

27. 

(7) 

20 

22 

.30 

the 

31 

the 

Ji 

in- 

33 

3‘i 

ill  y  -  y'cos'i {>  +  z'siivj) 

iil  z  =  -y'siri<j>  +  z'coniji  (9) 

Substituting  this  expression  into  Equation  (8)  and  collecting 
---  terms  of  integration  of  dy' ,  we  find  the  integral 

*l9_  ! 

50.  j 

5L  - f exp[  j(kvcos^---kzsin<t>)y'.  }dy.'.9;  =  ,2.ir5(kycos<!>  — kzsin<j>> — (-10)  i 

J  ■>  >■  ur  .  •*  I 
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EVE?i 


If  we  define 


we  can  write  the  Dirac-6  function  of  Equation  (10)  as  27i6(ky“ 
ky*).  Equation  (8)  can  now  be  written  as 


!^X^o-zXs_lZ.2ir"/.p0.Ckx.,.ky>6"Ck.yn;kyl)_X 


exp[  jCk^x+kyz’  sin$  +  k^z’cos^JJdkjjdky 


(12) 


The  integral  over  dy*  can  be  evaluated  using  the,  sifting _ proper.- 

ties  of  the^lol-fuiic’t ioril 4 > ■-  -v-  • _ ! _ 


AND  ADDRESS 


p(x,<J>,z)  =  1/.2tt  J*Po(kx»ky* )  ^ 


ft. 


exp  t  j  (  kxx+ky  *  z  *  sin<|)+kz  z '  coscji ) )  dkx 


(13) 


« 

Since  the  origin  of  the  coordinate  system  is  in  the  plane  of  the 
light  path  slice,  we  have  z'  =  z^  =  0.  Equation  (13)  now  becomes  j 

I 

I 

t 

P(x,<|>,z0)  =  l/2it  j*p0(kx,ky,)exp(jkxx)dkx  (14) 

1 

I 

J 

Equation  (11)  can  be  restated  as! 

(15a} 


ky*  -  (k2  j-  kx2 )sin<(>  . 


If  kx  is  small  compared  to  k,  then  the  following  approximation! 
holds.  i  .  t 


i 

ky'  =  ksin<£ 


(15bj 


1  7  •  i 

j-'y  U.'hen  this  approximation  is  substituted  into  Equation  (14), 
if  !p(x,<),z0)  becomes  p(x-,ky*,z0) 
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p(x-;ky* , z-0-y  =.  1/ 2tt  f  P0(.kx-,ky’  )exp(-jkxx)dkx - (16) 
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It  is  important  to  note  that  pCx.ky'^o)  lies  midway  between  _J> 
the  Fourier  domain  and  the  time-s'pace  domain.  If  we  take  the  in-  6 
verse  Fourier  transform  of  p(x,ky!’  ,zD)  with  respect  to  the  vari-  7 
able  kv‘ ,  we  obtain  !  8 

;  9 

./p.(x,kyt , z0 )'exp (  jky 1  y )dky *  « _ id 

_ L _  12" 


8-CKxj  ky: j  z0)^pj  iUixx±k  ♦  y )J  dkxdk 


The  term  o rv.jjt he> Iri^h tlha'nd.  s_ide. . _c an_bg_te_C.O gtJ.i zp.d __  as _ pCXvy.,.a0)..._|  l£ 

Thus  from  _a _ set _ of _ measurement  taken  at  a  given  elevation  (x  j  _17 

fixed)  and  .yaWink1.  angle';!  We^'can  japplv  a  one-dimensional  Fourier"]  1_8 
transform  _fip^.b.t.ain_.the..pr.essj4tLei_aiQ.ng_tl'ifi_.y^axi.s_f.oxJthat._YaXuc..]  JL2 

of  x.  The  .pMssy$e-J>yer  _a .  specified.  >ty_.pJ,a.ae_.can,.alao.  .be _ .Qbt.arJ  20 

ined  by  a  series  of  such  transforms  taken  at  equally  spaced  va-  J  21 
lues  of  „x .  • _ L _ 1  22 


If  we  take  the  Fourier  transform  of  p(x,ky’,z0>  with  respect 
to  the  variable  x,  we  obtain  ! 


1 

p(kx,ky’,z0)  x^ky*  ,z0)texP(-jkxx)Mkx 


L  which  is  the  pressure  field  in  the  Fourier  domain. 

?  ! 

—  t 

3.  Implementation  of  Method  B  using  DFFT  algorithms  requires 

!i.  the  approximation  that  ky’  does  rjot  depend  on  kx.  This  approxi- 
!>_  mation  is  valid  when  the  ultrasound  is  confined  to  a  narrow  beam 
!?_  as  with  sound  field  produced  by  most  medical  and  NDE  transducers. 

L  1  ! 

?  ! 

i~  To  illustrate  this  approximation  we  show  in  Figure  3  the 

nearly  rectangular  format  for  'data  obtained  from  Equation  (16) 

—  using  an  DFFT.  The  error  incurred  by  assuming  this  format  t:o  be 
L  rectangular  will  be  snail  if  most  of  the  radiated  energy  is  near 
--  the  origin  in  the  Fourier  domain.]  The  illustration  in  Figure  3 

—  jis  for  a  sound  field  produced  by  a  circular  transducer  of  radius 
y-  a  =  10X  where  X  is  the  acoustic  Wavelength.  Each  of  the  larger 

—  concentric  circles  has  associated  with  it  a  percentage  of  total 
?;■  radiated  energy  contained  within!  the  circle.  The  percentages 
i-  were  calculated  using  an  Airy  ipattern  to  approximate  the  sound 

field.  Figure  3  shows  the  format  to  be  essentially  rectangular 
~  for  98%  of  the  radiated  energy  in  this  case.  Figure  4  shows  the 
general  curve  for  the  fraction  of  total  energy  radiated  as  a 
J—  function  of  ka(sin<j>-)-  for- the.  Airy  , function.  - - - 


Distribution  of  radiated 
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4.0  DATA  acquisition: u sing:-' METHOD; 


I  I 

Equations  (15a)  and  (15b)  indicate  that  the  angle  $  should  j 
change  by  equal  increments  of  sin  $ ,  or  that’ 


sin<!>  =  +/-  not 


(19)  j 


where  n  =  0,1,2,...,N  with  2N+1  H  number  of  experimental _ points  1 

axidOgO=Ispecrfi'edTihcrement  f or  sin _ 


1 


1  t; 


The  maximum  value  of  k  observed  in  the  data  is 


From  Equation  (15b),  we  also  know 


(ky)max  ~  ,k(sin<j>max) 
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(21) 


Combining  Equations  (19),  (20)  arid  (21)  yields 


Akv  |=  ka 

+  I 


(22)  , 


The  increments  of  y  in  the  transformed  data  returned  by  the  ( 
DFFT  are  y  <=  mAy  where  m  =  0,1, 2 M  and  M  is  the  number  of  po-  ! 
ints  used  in  the  DFFT.  From  the  ^periodicity  of  the  DFFT,  we  find 


or 


MAkyAy  =  2it 

I 

I 

I 

I 

I 

I 

Ay  =  c/(Maf) 


(23) 


(24)  | 


where  c  is  the  sound  speed,  and  the  acoustic  frequency  f  *=  c/1.  | 

jThe  total  range  on  the  y-axis  is  then  j 

i 


MAy  <=  c/(ctf)  . 


(25)  | 
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5.0  EXPERIMENTAL  RESULTS 


zj 


Cook  and  Berlinghieri  have  experimentally  demonstrated  the  j 
validity  of  Method  A.^  They  generated  a  rectangular  array  of 
points  in  the  Fourier  domain  by  rotating  the  transducer  through 
the  proper  angles  and  sampling  the  linear  scan  at  proper  inter-  j 
vals.  They  then  used  a  two-dimensional  DFFT  to  reconstruct  the 
acoustic  field  in _ the _ plane_  of _ interrogation  and  in  nearby 


LHetho'd 3 ^wasTused  tq _ calculate. _ the _ pressure _ d is ,t xib.u.Uio ruJ 

over  a  transverse  plane  for  a  1;.25  cm.  diameter  PZT  transducer 

submersed  in._K3 .tfi x. _ Thaj: xans duce.x„v;a.s_.op.era  ting_.at_2 .2 -blhz._ (n_s_ 

9.3  A  ).  A^'i'iqgok'al U‘loffhq.1  *- 41  ;data-points..using-a_value-£)f.-_ij)-'v- 

sin(  1  degree).  .were_us.ed-J:.Q.jneas.uxe—the-4>xes.s.uxcL-in_a_.pIane_5 _ cm 

from  the  “  i.xapsdjucfitiiXf ace'i!-[^!fhi|s_data_.was. .used ..to.. cons tx uct_the  J 

pressure  f ield_o.v.ex_an_axea_3..9_x;_3..9_cm. _ .Tbe_results.are....shown.J 

|in  Figure  m  J-$fi§ie£f  ects„_of._  the._ap.pxon£imat±on  .made ,. in.. Equation  _. 
(15)  are  noticeable  in  the  reconstruction.  The  reconstructed- 
field. _ia  .slightly— elliptical, rather  than- circular- due -to- the _out--_! 
ward-displaaementLjblUti^Li^zw&ve-yectotrcomponerits  _  in-.the —  recon-  ! 
struction.  | 
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